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MODE  APPROXIMATIONS  FOR 
IMPULSIVELY  LOADED  RIGID-PLASTIC  STRUCTURES 

by 

1  2 
J.  B.  Martin  and  P.  S.  Symonds 


Summary: 

An  attempt  is  made  to  provide  a  rational  method  of  constructing 
one  degree  of  freedom  approximations  for  impulsively  loaded  metal 
structures  which  are  analysed  according  to  an  elementary  rigid-plastic 
theory.  The  approximation  follows  automatically  from  a  chosen  mode 
shape,  and  a  criterion  for  determining  good  mode  shapes  is  introduced. 


Assistant  Professor  of  Engineering,  Brown  University. 
Professor  of  Engineering,  Brown  University. 
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1.  INTRODUCTION 

The  problem  of  the  estimation  of  permanent  deformations  in  ductile  met 
structures  subjected  to  high  intensity  dynamic  loading  has  received  increasing 
attention  over  the  past  fifteen  years.  The  problem  is,  in  its  entirety,  one  of  gre 
difficulty.  Above  the  complexity  of  the  problem  of  elastic  structures  subjected  ti 
transient  loading,  additional  complicating  factors  which  are  of  importance  includ 
the  dissipation  of  energy  in  plastic  work,  elastic  unloading  from  plastic  states,  \ 
hardening,  the  dependence  of  yield  stress  upon  rate  of  strain,  geometry  change  t 
and  divers  other  non-linear  effects.  Very  few  solutions  to  dynamic  loading  prob 
have  recognized  all  these  factors.  More  commonly,  one  or  i  lore  of  the  factors  i 
assumed  to  dominate  the  behavior  of  the  structure,  and  all  others  are  neglected  ' 
approximated. 

Probably  the  most  widely  used  approximation^  is  the  replacement  of  the 
distributed  mass  of  the  structure  by  one  or  more  point  masses.  This  approxima 
has  been  used  in  conjunction  with  a  variety  of  idealizations  of  the  material  behav 
In  fairly  simple  structures  the  actual  distribution  of  mass,  elastic  stiff: 
and  yield  stress  can  be  considered.  Solutions  for  an  elastic,  perfectly  plastic 

^See,  for  example: 

N.  M.  Newmark,  "A  method  of  computation  for  structural  dynamics,”  Proc.  E 
Mech.  Div.,  A.S.C.E.,  (EM3),  p.  67,  1959. 

C.  H.  Norris,  et  al.  Structural  Design  for  Dynamic  Loads  ,  McGraw-Hill  Bool 
Co.,  New  York,  N.  Y.,  1959. 
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material  have  been  found  ,  but  are  difficult  if  the  load  magnitudes  are  much  larger 

than  those  causing  the  yield  stress  to  be  reached.  A  variety  of  dynamic  loading 

problems  has  been  solved  using  what  will  be  termed  in  this  paper  an  elementary 

rigid-plastic  theory.  A  representative  but  by  no  means  complete  list  of  papers  is 
5 

given.  The  elementary  rigid-plastic  theory  involves  the  following  idealizations: 
(i)  For  the  purpose  of  the  dynamic  loading  problem  under  consideration 

a  ductile  material  is  represented  by  a  rigid,  perfectly  plastic  constitu¬ 
tive  equation.  All  elastic  effects  are  in  consequence  excluded. 

(ii)  Geometry  changes  are  assumed  to  be  small,  and  the  yield  stress  is 
assumed  to  be  independent  of  the  rate  of  strain. 

4 

For  example: 

H.  H.  Bleich  and  M.  G.  Salvador!,  "Impulsive  motion  of  elasto-plastic  beams," 
Trans.  A.S.C.E.,  120,  p.  499,  1955. 

J.  A.  Seiler,  B.  A.  Cotter,  and  P.  S.  Symonds,  "Impulsive  loading  of  elastic- 
plastic  beams,"  J.  App.  Mech. ,  p.  515,  1956. 

3 

E.  H.  Lee  and  P.  S.  Symonds,  "Large  plastic  deformations  of  beams  under 
transverse  impact,"  J.  App.  Mech.  ,  19 ,  p.  308,  1952. 

M.  F.  Conroy,  "Plastic-rigid  analysis  of  long  beams  under  transverse  impact," 

J.  App.  Mech.  ,  J^,  p.  465,  1952. 

P.  S.  Symonds,  "Dynamic  load  characteristics  in  plastic  bending  of  beams," 

J.  App.  Mech.  ,  Z^,  p.  475,  1953. 

H.  G.  Hopkins  and  W.  Prager,  "On  the  dynamics  of  plastic  circular  plates," 
z.  angew.  Math,  u  Phys.,  p.  317,  1954. 

E.  W.  Parkes,  "The  permanent  deformation  of  a  cantilever  struck  transversely 
at  its  tip,"  Proc.  Roy.  Soc.,  A228  ,  p.  462,  1955. 

P.  G.  Hodge,  Jr.  ,  "Impact  pressure  loading  of  rigid-plastic  cylindrical  shells," 
J.  Mech.  Phys.  Sol.,  3,  p.  176.  1955. 
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The  use  of  the  elementary  rigid-plastic  theory  is  the  subject  of  this  paper. 
It  is  not  within  the  scope  o^  the  report,  however,  to  give  a  complete  appraisal  of 
the  utility  and  applicability  of  this  theory.  Itis  entirely  a  valid  method  of  approach, 
as  are  the  other  idealizations  mentioned  above,  provided  that  care  is  taken  to  define 
the  range  where  idealization  and  reality  have  something  in  common,  and  provided 
that  experimental  evidence  confirms  the  theoretical  predictions  in  this  range.  The 
elementary  rigid-plastic  theory,  when  applied  to  structures  subjected  to  short 
duration,  high  intensity  loading,  has  a  range  of  validity  bounded  at  one  end  by  the 
requirement  that  the  deformations  should  not  be  so  large  that  geometry  change 
effects  are  significant,  and  at  the  other  end  by  the  requirement  that  the  energy  of  thi 
disturbance  should  be  large  compared  to  the  energy  which  could  be  stored  elasticall 
in  the  structure  in  order  to  justify  the  exclusion  of  elastic  effects.  These  require¬ 
ments  are  in  some  senses  contradictory  since  the  large  disturbance  needed  to  meet 
the  energy  requirement  will  tend  to  produce  large  deformations.  Thus  the  extent 
of  the  possible  range  of  validity  of  the  elementary  rigid-plastic  theory  must  depend 
on  the  configuration  and  flexibility  of  the  structure,  and  in  some  cases  may  not  exis 
at  all. 

Experiments  show  that  the  rigid-plastic  theory  almost  always  requires 
corrections  for  the  dependence  of  yield  stress  on  strain  rate,  strain  hardening  and 
geometry  changes.  When  these  corrections  are  made,  remarkably  good  agreement 
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has  been  achieved  ,  showing  that  the  neglect  of  elastic  effects  is  permissible  when 

the  total  energy  dissipated  is  much  greater  than  maximum  energy  which  could  be 

stored  elastically,  '^n  appropriate  circumstances  it  is  possible  to  make  the 

corrections  for  strain  rate  sensitivity  and  for  strain  hardening  by  simply  multiplying 

7 

the  static  yield  stress  by  a  constant  factor  ;  however,  the  "appropriate  circumstance 
are  not  always  obvious,  and  this  method  should  be  used  with  caution  for  a  highly  rate 
sensitive  material  such  as  mild  steel. 

In  the  authors'  opinion  the  importance  of  the  elementary  rigid-plastic  theory 
lies  in  its  ability  to  provide,  quickly  and  simply,  an  estimate  of  major  deformations 
due  to  very  large  dynamic  loads.  Such  an  estimate  provides  a  convenient  basis  for 
more  refined  analyses  and  calculations  to  include  effects  of  strain  rate  sensitivity, 
finite  deflections,  and  other  effects  when  necessary. 


S.  R.  Bodner  and  P.  S.  Symonds,  "Ejq>erimental  and  theoretical  investigation  of 

the  plastic  deformation  of  cantilever  beams  subjected  to  impulsive  loading,  " 
J.  App.  Mech.  ,  29,  p.  719,  1962. 

E.  H.  Lee  and  S.  J.  Tupper,  "Analysis  of  plastic  deformation  in  a  steel  cylinder 
striking  a  rigid  target,  "  J.  App.  Mech.  ,  21,  p.  63,  1954. 

T.  J.  Mentel,  "The  plastic  deformation  due  to  impact  at  a  cantilever  beam  with 

an  attached  tip  mass,  "  J.  A.  M.  ,  p.  515,  1958. 

T.  C.  T.  Ting  and  P.  S.  Symonds,  "Impact  of  a  cantilever  beam  witli  strain  rate 
sensitivity,"  Proc.  Fourth  U.  S.  Nat.  Congr.  App.  Mech.,  A.  S.  M.  E.  , 
p.  1153,  1962. 

A.  J.  Frick  and  J.  B.  Martin,  "The  plastic  deformation  of  a  bent  cantilever," 

Tech.  Kept.  NSF  GPU  15/16  from  Brown  Un  iversity  to  the  National 
Science  Foundation,  April  1965. 


W.  Parkes,  "The  permanent  deformation  of  a  cantilever  struck  transveisely  at 
its  tip,"  Proc.  Roy.  Soc. ,  A  228,  p.  462,  1955. 
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This  report  will  be  limited  to  the  discussion  of  one  dimensional  coni 
i.e,  ,  structures  made  up  of  bars  and  rods,  curved  or  straight,  whose  cross 
dimensions  are  small  compared  with  their  length.  Generalized  stresses  (bei 
moments,  axial  force,  etc.)  and  generalized  strains  (curvatures,  axial  strai 
will  be  used  in  the  analysis.  The  following  description  of  the  plastic  flow  ru] 

g 

rigid-perfectly  plastic  material  follows  Prager. 

Let  the  generalized  stresses  acting  at  a  section  be 

Qj  (j  =  1. ...  .  n) 

and  let  the  associated  generalized  strain  rates  be 
q.  (j  s  1, ....  n). 

J 

The  dimension  of  a  particular  component  of  generalized  strain  is  such  that  tl 
stress-strain  product  has  the  dimension  of  work  per  unit  length.  Tuc  stress 
relation  is  written  in  terms  of  a  yield  function  '1*{Q.)  which  must  be  convex  a 
which  must  contain  the  origin.  The  generalized  strain  rates  are  given  in  ter 
the  generalized  stresses  as  follows: 

q.  =  X<  *(0.)  >-5^  (1) 


where 


<  *(Q^)  >=  1 


when  <!'  (Q .)  =  0 


W.  Prager,  An  Introduction  to  Plasticity.  Addison  Wesley  Press,  (Readinj 
Massachusetts),  1959. 
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<  *(0.)  >=  0  when  *(Q.)  <0 

X  >  0  but  otherwise  unspecified. 

Stresses  such  that 

♦  (Q.)  >0 

are  not  admitted. 

The  geometric  interpretation  of  equation  (1)  is  now  well  known.  If  the  r 
components  of  the  stress  are  plotted  as  coordinates  in  an  n  dimensional  space 
a  stress  point  is  defined.  Similarly  the  convex  function  <l>  =  0  may  also  be  plott' 
as  a  closed  surface  in  the  stress  space.  If  the  stress  point  lies  inside  the  surfac 
♦  =  0,  the  strain  rate  is  zero.  If  the  stress  point  lies  on  the  surface,  the  magnit 
of  the  strain  rate  vector  is  not  specified.  However,  if  q.  is  plotted  in  the  stress 
space  with  the  stress  point  as  origin,  the  strain  rate  vector  is  required  to  have  t 
direction  of  the  exterior  normal  to  the  surface  0  at  that  point.  A  two  dimens 
representation  of  tliis  interpretation  is  given  in  Fig.  1. 

♦ 

This  geometrical  approach  demonstrates  that  if  Q.  is  an/  other  admis 
state  of  stress  (Fig.  1),  then 

(Q  -  Q  *)q  >0  .  (2) 

J  J  J 

♦ 

Further,  if  itself  lies  on  the  yield  surface  and  is  associated  with  strain  rate 

.  * 

q^  (shown  in  Fig.  1),  it  is  clear  that  by  a  second  application  of  the  concept 
deomonstrated  in  equation  (2), 

(Q.  -  Q.*)  (q.  -  q  *)  >0 

J  J  J  J  - 


(3) 
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These  inequalities  have  been  discussed  in  a  more  general  context  by  Hill  and 
Drucker. 

Inequality  (3)  may  be  used  to  show  that  the  velocity  history  at  any  point 

in  a  structure  following  the  application  of  high  intensity  time  dependent  loading  is 

unique,  provided  that  changes  in  geometry  may  be  ignored  and  that  ce-  tain  continuity 

requirements  are  satisfied.  This  result  has  been  discussed  by  Martin^  \  and  is 

12 

based  upon  earlier  work  by  Drucker  which  appears  to  have  passed  unnoticed  by 
most  workers  in  this  field. 

Suppose  that  a  given  structure  is  subjected  to  the  following  loading  and 

I 

boundary  conditions:  On  length  S^of  the  body  time  dependent  tractions  T. 

(i  =  1,  2,  3)  are  specified,  and  on  the  remainder  of  the  structure  S  time  dependent 

u 

t 

velocities  u^  are  prescribed.  S^^and  S^,  which  together  comprise  the  whole 
structure,  may  themselves  be  time  dependent.  Further,  let  the  velocities  at  time 
t  s  0  be  given  by  v^.  These  quantities  define  the  problem.  Assume  now  that  two 
solutions  can  be  found.  First,  velocities  u.  associated  with  accelerations  u.  , 


R.  Hill,  "A  variational  principle  of  maximum  plastic  work  in  classical  plasticity,  " 
Quart.  J.  Mech.  App.  Math. ,  1,  p.  18,  1948. 

10 

D.  C.  Drucker,  "A  more  fundamental  approach  to  stres s -strain  relations," 

Proc.  1st  U.  S.  Congr.  App.  Mech.,  A.S.M.E.,  p.  487,  1951. 

B.  Martin,  "A  note  on  the  uniqueness  of  solutions  for  dynamically  loaded 

rigid-plastic  and  rigid-visco-plastic  continua,"  Tech.  Rept.  NSF  GPU  15/12 
from  Brown  University  to  the  National  Science  Foundation,  September  1964. 

12 

D.  C.  Dricker,  "A  definition  of  a  stable  inelastic  material,"  J.  App.  Mech.  ,  26, 
p.  101,  1959. 
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resses  Q  ,  strain  rates  q..  Secondly,  velocities  ii.  ,  accelerations  u,  ,  stresses 

J  J  1  1 

,  strain  rates  .  Both  solutions  satisfy  all  the  field  equations. 

Consider  the  two  solutions  at  some  time  t  >t  .  Since  the  tractions,  inertia 

o 

rces  and  stresses  of  each  solution  are  in  internal  and  external  equilibrium,  tlie 

'ferences  in  these  quantities  are  in  equilibrium.  Similarly,  the  differences  in  the 

rain  rates  and  velocities  are  compatible.  Further,  the  difference  in  the  tractions 

nishes  on  S_i  and  the  difference  in  displacements  vanishes  on  S  .  In  consequence, 
1  u 

>m  the  principle  of  virtual  velocities, 


-  u.  )  (u.  -  u.  )  dS  = 
11  1 


(Q.  -  Q,*)  (q.  -  q."')  dS 


(4) 


ere  m  is  the  mass  per  unit  length  of  the  structure.  Fron,  the  right-hand  sidi’ 
non-negative.  Hence,  provided  that  the  velocity  at  each  point  is  a  continuous 
iction  of  time,  we  may  write 


-  u.  )  dS  <  0 
1 


(5) 


is  seen  that  a  non-negative  quantity  involving  the  ve  locities  u.,  ii.  must  decn  as*- 

jje 

th  time.  At  time  t  =  0,  however,  ii.  =  u.  =  v.,  and  the  non-negative  quantity  is 

. 

tially  zero.  It  follow  s  that  u.  =  u  for  all  t  >0. 

11 

The  restriction  imposed  by  the  requirement  that  the  velocity  at  each  point 
ould  be  a  continuous  function  of  time  is  a  serious  one  in  general.  It  may  readily 
shown  that  when  discontinuities  are  permitted  the  solutions  to  certain  problems  are 
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not  unique  (Ting  ).  In  practice,  however,  a  large  class  of  problems  where 
discontinuities  do  not  occur  is  of  interest.  This  is  the  class  of  structural  problems 
where  shear  and  axial  strains  are  stipulated  zero  (bending  and  torsion  strains  are 
permitted).  The  velocity  at  each  point  will  be  a  continuous  function  of  time  in  these 
cases,  and  hence  uniqueness  follows.  This  report  will  be  restricted  to  problems 
falling  into  this  class. 

Further,  the  report  will  be  concerned  only  with  impulsive  loading  problems 
These  problems  may  be  characterized  in  exactly  the  same  way  as  the  general 
problem  given  above  with  the  following  restrictions: 

(i)  The  tractions  applied  to  will  be  taken  to  be  zero. 

(ii)  The  velocities  prescribed  on  S  will  be  taken  to  be  zero. 

u 

(iii)  and  will  be  assumed  to  be  time  independent. 

The  impulsive  loading  problem  is  thus  essentially  one  in  whicl'.  initial  velocities  are 
prescribed  over  the  whole  structure  at  time  t  =  0;  tliereaftcr,  no  external  forces  do 
work  on  the  structure. 

No  direct  technique  exists  for  determining  the  response  of  the  structure  to 
this  (or  any  other  more  complex)  form  of  loading.  A  solution  must  be  sought  on  a 
trial  and  error  basis.  If  a  solution  can  be  found  which  satisfies  all  the  conditions 
(equilibrium,  compatibility,  yield  and  the  stress  -  strain  relation),  the  uniqueness 
result  assures  that  this  is  the  only  solution  for  the  velocity  as  a  function  of  space 
and  time. 


13 


T.  C.  T.  Ting.  Private  communication,  March  1965. 
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The  response  of  a  one-dimensional  rigid-plastic  structure  to  impulsive 
>ading  is  characterized  by  two  distinct  phases  of  behavior.  In  the  first  phase, 
ravelling  "hinges"  (discrete  points  at  which  plastic  deformation  occurs)  associated 
ith  discontinuities  in  the  acceleration  field  are  found.  In  the  second  phase,  deforma - 
on  occurs  without  change  of  shape  of  the  velocity  field.  The  equations  of  the  second 
base  may  be  written  in  terms  of  separate  functions  of  time  and  space  parameters, 
nd  hence  are  easily  formulated. 

14 

In  considering  the  elementary  rigid-plastic  solutions,  Symonds  emphasized 

lat  in  certain  cases  very  simple  closed  form  solutions  of  impulsive  loading  problems 

an  be  obtained,  for  example  by  momentum  conservation  equations  in  finite  form, 

15 

xamplcs  were  g.  ’en  where  slightly  modified  problems  required  numerical  irtegra- 

.on  of  the  equations.  Further  examples  have  been  given  in  the  literature.  Thus  the 

roblem  of  a  mass  striking  a  fixed  end  beam  was  solved  in  closed  form  by  Parkes 

ut  for  the  pin  ended  case  the  equations  for  the  first  phase  have  to  be  solved 
17 

umerically  (Ezra  ). 


P.  S,  Symonds,  'Large  plastic  deformations  of  beams  under  blast  loading,  ” 

Proc.  Second  U.  S.  Nat.  Congr,  App,  Mcch.  ,  A.S.M.E.,  p.  505,  1954. 

) 

P.  S.  Symonds,  "Simple  solutions  of  impulsive  loading  and  impact  problems  of 
plastic  beams  and  plates,"  Tech.  Kept,  UERD-3  from  Brown  University 
to  the  Norfolk  Naval  Shipyard  under  Contract  Ni895-1756  A,  April  1955. 

*E.  W.  Parkes,  "The  permanent  deformation  of  an  encastre  beam  struck 

trans-'ersely  at  any  point  in  its  span,"  Proc.  Inst.  Civ.  Engrs.  (London), 
10,  p.  266,  1958. 

A,  A.  Ezra,  "The  plastic  response  of  a  simply  supported  beam  to  an  impact 

load  at  the  center,"  Proc.  3rd  U.  S.  Nat.  Congr.  App.  Mech.  ,  A.S.M.E., 
p.  513,  1958. 
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A  further  example  is  given  in  Fig.  3.  This  is  a  fixed  end  beam  subjected  to 
impulsive  loading  so  that  the  velocity  is  in  the  region  of  length  b.  The  complete 
solution  can  be  easily  obtained  in  closed  form  for  any  ratio  b/i.  However,  if  the 
ends  are  pinned  rather  than  clamped  a  quite  unpleasant  numerical  solution  is  require 
in  order  to  obtain  the  three  quantities  required  to  define  the  motion  in  the  first  phase 
This  example  will  be  discussed  later  in  the  paper. 

Although  the  solution  by  numerical  integration  of  a  problem  of  elementary 
rigid-plastic  theory  is  almost  certain  to  be  much  less  difficult  than  the  solution  of 
the  same  problem  according  to  elastic -plastic  theory,  it  may  be  tedious  enough  so 
that,  for  purposes  of  first  estimates  at  least,  one  would  welcome  a  simple  approxim 
method.  Such  methods  have  frequently  been  used.  The  most  common  method  is 
equivalent  to  the  replacement  of  the  structure  by  a  mass -spring  system  of  one  degre 
of  freedom.  Matching  the  model  to  the  actual  problem  has  been  done  in  various  ad 
hoc  ways.  No  attempt  seems  to  have  been  made  to  generalize  or  systematize  the 
various  approximate  methods. 

This  report  will  discuss  a  technique  of  approximation  which  is  closely 
related  to  the  matching  of  one  degree  of  freedom  models.  Using  an  argument  akin 
to  that  used  to  establish  uniqueness,  it  will  be  shown  that  the  solution  to  one  impulsi’ 
loading  problem  may  be  approximated  by  the  solution  of  a  problem  involving  the  sam 
structure  but  different  initial  velocities. 

To  make  use  of  this  property  we  shall  first  give,  for  illustrative  purposes, 
the  solutions  to  two  simple  problems,  and  then  show  how  approximations  may  be 
constructed  for  these  and  other  problems.  Finally,  we  shall  discuss  the  relation 
between  the  two  approaches  to  approximation,  showing  that  the  concepts  of  the  previ< 


I 
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paragraph  may  be  used  to  provide  a  rational  method  of  matching  one  degree  of 
freedom  models  to  the  structure. 
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2.  SOME  EXAMPLES  OF  SIMPLE  SOLUTIONS 

In  this  section  we  shall  present  two  problems  which  can  be  very  simply 
solved  with  a  rigid-plastic  theory.  These  solutions  are  given  primarily  to  illustr 
problems  where  the  solution  for  the  first  phase,  when  travelling  hinges  occur,  ca 
written  in  closed  form 

(i)  The  first  problem  considered  will  be  the  uniform  contilever  subjected  to 

transverse  impulse  on  a  point  mass  attached  at  the  tip.  The  solution  to  this  prob 

7 

was  obtained  by  Parkes  .  The  beam,  with  the  relevant  physical  parameters,  is 

shown  in  Fig.  4.  At  time  t  =  0,  the  tip  mass  acquires  a  velocity  v.  Shear  strai 

are  neglected  in  the  analysis,  and  hence  the  bending  moment  is  constrained  to  be 

between  +M  and  -M  for  all  shear  force  values, 
o  o 

In  accordance  with  the  trial  technique  for  solution,  a  response  is  postula 

and  then  checked.  An  assumed  velocity  field  at  time  t  >  0  is  shown  in  Fig.  5(a). 

The  section  AB  of  the  beam  rotates  about  B,  and  thus  the  bending  moment  at  B  m 

be  zero.  3C  remains  stationary,  and  thus  for  equilibrium  of  BC  [Fig.  5(b)]  the  s 

force  in  BC  must  be  zero  and  the  bending  moment  must  be  M  .  Certainly  this 

o 

distribution  of  velocity  and  strain  rate  is  kinematically  admissible,  and  the  stres 
strain  relation  has  been  satisfied. 

We  give  the  solution  obtained  by  Parkes^  by  writing  conservation  of 
momentum  equations  for  the  beam.  Since  no  external  forces  act  on  the  beam  the 
linear  momentum  of  the  beam  cannot  be  changed:  hence  at  time  t 

Mxz 


Gz  + 


2 


=  Gv 


(6) 
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Angular  momentum  about  a  line  through  A  perpendicxilar  to  both  the  beam  center  line 
and  the  direction  of  motion  will  be  considered.  The  initial  angular  momentum  about 
this  line  is  zero;  hence  at  time  t 


,  mx  . .  X  . .  ^ 

These  equations  give  z  and  t  in  terms  of  the  parameter  x.  Introducing  the  dimension 

0 

less  parameters 


Y  s  mf/2G  ,  C  =  3c/i 

(6)  and  (7)  become 

(-£— )  - i - 

'  V  '  1  +  Y? 


(8a) 


M  t  .  J. 

o  1  rC 

Giv  "3  (1  +  y5) 


(8b) 


5  increases  monotonically  to  give  increasing  t,  and  the  solution  remains  valid  until 
^  =  1  or  the  hinge  reaches  C.  Thereafter  the  beam  rotates  as  a  rigid  body  about  C, 
and  the  behavior  enters  the  second  phase.  The  shear  force  at  C  is  no  longer  zero: 
however,  angular  momentum  about  C  gives  a  relation  between  z  and  t  direc  tly: 


M  t  =  -X- 
o  3 


Giv 
(1  +Y) 


+  Giv 


(1 

(1  +Y) 


-  Giz  (1  +-J-  Y) 


(9) 


Rearranging 


M  t 
o 


1  -  (-f-)(l  + 


3 


Y  ) 


Giv 


(10) 
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This  relation  applies  until  z  =  0. 

The  equations  derived  here  can  be  shown  to  lead  to  a  distribution  of  bendin 
moments  which  at  all  times  satisfies  the  yield  condition,  and  hence  this  is  the 
solution  to  the  problem.  A  closed  form  solution  for  the  final  displacement  6  of  the 
tip  mass  may  be  obtained  by  integrating  the  tip  velocity  equations: 


M  6 


o 


(H) 


(ii)  As  a  second  example,  consider  a  uniform  beam  subjected  to  impulsive 

loading  over  a  lengtli  b  symmetric  about  the  midpoint,  so  that  the  initial  velocity  is 
v^  on  this  section  and  zero  in  the  remaining  length  {t  -  b)  [Fig.  6(a)].  The  yield 
condition  is  taken  to  be  the  same  as  that  in  the  previous  example.  For  b  =  i  the 
solution  was  given  by  Symonds  for  clamped  end  and  pin  end  conditions  ’  .  In 

that  case  it  is  not  difficult  to  obtain  solutions  in  closed  form.  In  practice,  if  the 
ends  are  constrained  against  axial  motion,  the  effects  of  the  constraints  tend  to 

dominate  when  the  deflection  exceeds  magintudes  of  the  order  of  the  beam  depth. 

18 

These  effects  have  been  considered  ,  but  they  lie  beyond  the  scope  of  the  present 
discussion. 

For  b/i  <  1,  if  the  ends  are  built  in  a  closed  form  solution  can  be  simply 
derived.  Fig.  6(c)  indicates  the  trial  velocity  configuration  for  an  instant  soon 
after*  motion  begins.  The  sections  AP  and  QC  have  zero  shear  force  and  bending 
moment  -M  and  M  ,  respectively.  The  plastic  boundary  points  are  at 


P.  S.  Symonds  and  T.  J.  Mentel,  "Impulsive  loadings  of  plastic  beams  with  axial 
constraints ,"  J.  Mech.  Phys.  Sol.,  6,  p.  186,  1958. 


X  =  (t)i 


and 

X  —  > 

and  these  define  the  motion  in  this  stage.  They  can  be  found  by  writing  the  two 
equations  of  momentum  conservation,  namely  of  linear  momentum  and  angular 
momentum  about  any  point  in  AP  or  QC.  Omitting  details,  the  following  results 
are  found: 


1 


2 


12  M  t 
o 


mi  V 


o 


(12) 


This  stage  ends  when  the  hinge  at  P  reaches  the  support  point  or  that  at  Q  reaches 
the  midpoint.  If  b  <  i/2,  the  subsequent  stage  involves  the  pattern  shown  in  Fig.  6(d), 
the  unknown  quantities  defining  the  motion  now  being  the  midpoint  velocity  v(t)  and 
Cj(t).  Momentum  equations  enable  us  to  find  these  as  follows: 


12  M  t 

£  ,  _S_  +  - 2. 

^1  4i  mbiv 

o 


4  V 


=  3  + 


48  M  t 
o 

mb^v 


o 


These  equations  hold  for  times  t  such  that 


1  < 


48  M  t 
o 

mb^v 


(13a) 


(13b) 


(14) 
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The  latter  time  corresponds  to  the  instant  when  the  hinge  at  P  reaches  end  A  of  th 
beam.  The  two  stages  treated  above  comprise  the  first  phase  of  the  motion. 

In  the  second  phase  there  are  stationary  plastic  hinges  at  A  and  C.  Eac^ 
half  beam  rotates  as  a  rigid  bar.  The  midpoint  velocity  v  is  found,  by  writing  an 
angular  momentum  equation,  to  be  given  by 


48  M  t 
o 


mb  V 

o 


The  beam  comes  to  rest  when  t  =  t 


f 


,  where 


48  M  t . 
o  f 

mb^v 

o 


3 


(15) 


(16) 


The  deflections  in  the  various  stages  are  easily  found  by  integrating  the  velocities 
written  above.  The  result  for  the  final  midpoint  displacement  is 


u. 


mb^v 


12  M 

o 


+ 


+  1).  b  < 


(17) 


The  three  terms  show  respectively  the  contributions  from  the  three  stages  of  mot: 

Throughout  the  above  analysis  the  yield  condition  is  everywhere  satisfied 
Since  the  equations  of  dynamics,  the  boundary  conditions  and  the  initial  velocity 
distribution  are  all  satisfied,  the  result  is  the  exact  solution  (according  to  the 
elementary  rigid-plastic  theory).  The  above  results  apply  for 


2b/i  <  1  . 
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A  similar  analysis  for 

2b/i  >  1 


gives  the  result 


mb^v 


48  M 

o 


(18) 


When  the  ends  of  the  beam  are  pinned  rather  than  built  in,  the  solution  for 

general  h/l  <  1  cannot  be  obtained  in  the  simple  manner  used  above.  The  velocity 

pattern  and  free  body  diagram  for  the  first  stage  are  shown  in  Figs.  7(b)  and  7(d). 

These  enable  the  initial  velocity  conditions  and  all  other  requirements  to  be  satisfied. 

There  are  now  four  unknowns;  in  addition  to  and  4^  velocity  v^  at  the  hinge 

point  P  and  the  reaction  force  must  be  found.  The  four  equations  in  these 

unknowns  can  be  written  in  various  ways  and  will  be  omitted  here.  Velocities  are 

continuous  at  the  hinge  points  P  and  Q,  but  the  accelerations  are  not.  The  numerical 

integrations  required  to  obtain  the  complete  solution  are  tedious.  They  have  been 

1 9 

carried  through  for  h/l  =  0.2  and  h/l  =  0.7. 

The  two  examples  given  show  just  how  easily  some  solutions  in  the  elemental 
rigid-plastic  theory  may  be  obtained.  The  case  discussed  in  the  previous  paragraph 
shows,  on  the  other  hand,  how  a  small  change  in  initial  or  boundary  conditions  may 
enormously  complicate  the  analysis.  Thv'  need  for  a  simpler  way  of  obtaining  an 
approximation  to  the  complete  solution  in  such  cases  is  apparent. 

19 

G.  Gangopadliyay .  M.  Sc.  Thesis,  Brown  University,  1964.  Further  calculations 
on  this  problem  have  been  carried  out  by  K.  Vashi. 
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3.  ONE  DEGREE  OF  FREEDOM  APPROXIMATIONS 

A  one  degree  of  freedom  model  appropriate  to  a  rigid-plastic  material  is 
shown  in  Fig.  8(a).  The  spring  force-spring  displacement  characteristic  is  as 
shown  in  Fig.  8(b);  essentially  motion  of  the  mass  is  resisted  by  a  constant  force  S 
If  the  spring  is  initially  undeformed  and  has  initial  velocity  u^,  it  is  clear  that  the 
velocity  u  at  tim.e  t  is  given  by 


Gix  =  Gu  -  St 
o 


for  u  >  0 


or 


St 

Gu 

o 


for 


(19) 


Thus  the  velocity  of  the  mass  decreases  with  time,  coming  to  rest  in  time  t^  =  Gu 
The  final  displa  ement  is 


u.  =  Gu  ^/S. 
f  c  o 


In  order  to  apply  these  results  to,  for  example,  a  beam  problem,  suitablt 
values  must  be  found  for  the  quantities  u,  S  and  G.  These  quantities  may  be  refer 
to  as  the  equivalent  velocity,  equivalent  resisting  force  and  equivalent  mass 
respectively. 

Consider,  for  example,  the  beam  shown  in  Fig.  9(a).  This  case  was 
discussed  in  detail  in  the  previous  section.  A  one  degree  of  freedom  approximatio 
will  be  given  for  this  example  in  order  to  demonstrate  a  method  of  finding  the 
equivalent  quantities. 
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In  the  case  of  the  spring  (Fig.  8)  the  limiting  value  of  the  resisting  force 
may  be  found  by  statically  applying  an  external  force  to  the  spring  and  finding  the 
value  of  the  external  force  required  to  initiate  flow.  This  external  force  and  the 
resisting  force  will  hav:;  equal  magnitudes.  This  simple  idea  provides  on('  vr.  tans 
finding  an  equivalent  force:  Apply  static  loading  to  the  beam  and  take  the  limiting 
value  of  the  load  parameter  as  the  equivalent  resisting  force. 

In  this  case  choose  the  loading  pattern  shown  in  Fig.  9(b).  Only  half  the 
beam  is  shown  because  of  the  symmetry  of  the  system.  Statics  shows  that  the 
limiting  value  of  P  is  given  by 


8  M 

P  - - 2_ 

c  2  i  -  b 


(20) 


Steady  flow  will  occur  in  the  beam  when  P  reaches  P  .  The  flow  field  associated 

c 

with  P^  suggests  itself  as  a  reasonable  means  of  choosing  the  equivalent  velocity; 
when  impulsive  loading  is  applied,  it  will  be  assumed  that  the  velocity  field  has  th 
same  shape  as  the  velocity  field  si  own  in  Fig.  9(c).  (Notice  that  the  velocity  fu  ld 
for  steady  flow  under  loads  P^  is  not  unique.  The  field  given  is  a  possible  con¬ 
figuration).  It  remains,  therefore,  to  determine  the  amplitude  of  the  velocity  fielc- 
as  a  function  of  time. 

One  method  of  determining  the  acceleration  is  to  write  the  angular 

.  Jje  i{c 

acceleration  equation  for  half  the  beam  about  the  hinge  support.  If  u  and  u  are 
the  velocity  and  acceleration  of  the  center  of  the  beam. 


mi 


u 


■)(-^) 


=  -2M 

o 


{ 


24 


(21) 
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Integrating  and  substituting  from  equation  (20) 


.*  *  6  (2f  -  b)  „  ^ 

U  s  V  -  . . . .  P  t 


where  v  is  the  initial  value  of  u  .  The  equivalent  mass,  in  analogy  with  equatio 


( 19),  is  given  by 


6  (2f  -  b) 


The  maximum  central  displacement  is  given  by 


.2  .  2 
mi  u 

o 

“f  ■  12  P  (2f  -  b) 
c 


It  is  now  necessary  to  ascribe  a  value  to  v  .  This  may  be  done  in  a  vari 
of  ways;  it  is  possible  to  match  the  initial  energy  of  the  actual  problem  and  the  on 
degree  of  freedorri  system,  or  to  match  either  the  linear  or  angular  momentum  of 
two  systems.  Computations  based  on  these  three  criteria  show  that  the  best  appr< 
imation  is  obtained  when  angular  momentum  of  half  the  beam  about  the  supporting 
hinge  is  balanced.  Thus  we  put 


mbv  . 


3b  (2i  -  b) 


V 


2 


(25) 


Substituting  into  (24)  and  rearranging  in  the  form  given  in  equation  (18),  the 
approximate  central  displacement  is  given  by 


(26) 


The  derivation  of  this  expression  involves  many  separate  decisions  which 
appear  natural  but  are  also  arbitrary.  If  the  method  is  to  be  applied  to  find  approx¬ 
imate  solutions  of  problems  which  have  not  been  fully  solved  it  is  desirable  that  the 
application  of  the  method  should  require  fewer  arbitrary  decisions.  A  more 
sophisticated  approach  using  the  principle  of  virtual  velocities  will  ^limunate  most 
of  these  decisions. 

This  method  is  based  on  the  premise  that  the  choice  of  a  mode  shape  is 
the  most  important  part  of  the  approximation  process.  Suppose  that,  tor  a  given 
structure,  the  solution  is  assumed  to  be  of  the  form 

u*  (S,  t)  =  S,  (S)  T  (t)  (27) 

where  S.  is  a  vector  valued  function  of  the  space  variable  S  alone,  and  T  is  a  scalar 

I 

function  of  time  alone.  Solutions  of  this  form  will  be  referred  to  as  mode  solutions. 

Some  freedom  is  permitted  in  arranging  the  magnitude  and  dimensions  of  the 

functions  S  and  T.  For  convenience  we  shall  assume  that  T  is  dimensionless  and 
i 

that  it  has  the  value  unity  when  t  =  0.  S^  will  thus  have  the  dimensions  of  velocity; 
it  will  in  fact  be  the  initial  velocity  distribution.  For  simplicity  in  solution  we  shall 
write  (27)  as 

u.*  (S.  t)  =  V.*  (S)  T  (t) 


(28) 
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If  the  solution  to  an  impulsive  loading  problem  is  to  be  assumed  to  be  of 

♦ 

the  form  of  (28)  it  is  clearly  desirable  that  (S)  should  satisfy  the  kinematic 
boundary  conditions. 

We  shall  show  first  that  T(t)  must  be  a  linear  function  of  time  if  the  intern 
energy  dissipation  rate  associated  with  the  assumed  velocity  field  is  to  be  equal  to 
the  rate  of  change  of  kinetic  energy.  The  rate  of  dissipa*^ion  of  the  kinetic  energy 
K  is  given  by 


m  .  *  .  * 
—T—  u.  u.  dS 
2  1  1 


=  L-''!  ''i 


dT  I  «  « 

-r-—  \  m  V.  V.  dS 

Js  ■  ‘ 


The  strains  q.  may  be  derived  from  the  velocities  u.  .  Let  the  initial  values  of 
J  1 

♦  r 

these  strains,  derived  from  v^  ,  be  q^.  These  strain  quantities  will  be  related  in 
the  same  way  as  the  velocities 


•  *  •  rr. 

q.  =  q.  T 
J  J 


When  the  strain  component  q^  is  plotted  as  a  vector  in  the  stress  space  (as  in  Fig. 
.  ’S' 

it  is  clear  that  q.  cnanges  only  in  magnitude  as  t*me  progresses.  It  is,  however, 

J 

.  *  ♦ 
the  direction  of  q^  which  determines  the  associated  stress  vector  Q.  .  (Although 

Q.  may  not  always  be  uniquely  determined,  the  dissipation  rate  is  unique  for  a 


given  q.  .)  Thus  Q.  does  not  change  with  time,  and  the  dissipation  rate  D  is  given 
J  J 

by 

«  *  .  ♦  «  r 

D  r  Q.  q.  =  Q.  q.  T  =  D  T  (31) 

J  J  J  J 

-  *  r 

where  D  =  Q.  q.  is  the  initial  dissipation  rate. 

J  J 

We  now  equate  the  rate  of  change  of  kinetic  energy  and  the  rate  of  change 
of  internal  energy,  given  by  the  total  energy  dissipation  rate  with  a  negative  sign. 


m  V.  V  dS  = 
1  1 


D  dS 


dT 

dt 


D  dS 


«  * 

m  V  V.  dS 
1  1 


X 


(32) 


X  is  a  constant,  defined  by  (32),  which  can  be  readily  de  ermined  from  the  initial 
velocity  distribution.  Solving  equation  (32)  for  T,  with  the  requirement  that  T  =:  1 
when  t  =  0,  we  obtain 


T  =  (1  -  Xt) 


(33) 


The  time  t^  which  elapses  before  motion  in  the  assumed  mode  ceases  is  given 
simply  by 


t 


* 

f 


(34) 


It  is  still  necessary  to  determine  the  initial  velocity  field  v.  ,  since  only 


the  mode  shape  and  not  the  initial  amplitude  has  been  assumed.  Suppose  that  a  mode 
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shape  is  assumed. 


Then 


V  *  s  o$.  (35) 

i  1 

and  it  is  necessary  to  choose  the  amplitude  a.  The  actual  initial  velocity  distri 

will  be  taken  to  be  v. . 

1 

A  rational  method  of  determining  a  may  be  developed  by  assuming  that 
one  degree  of  freedom  model  is  subjected, ’at  time  t  =  0,'to  an  impulse  distribul 

mv.,  where  m  is  the  mass  per  unit  length  of  the  structure.  The  velocity  with  w 

* 

the  model  begins  to  move  may  be  then  taken  to  be  v.  . 

This  problem  may  be  handled  by  means  of  the  principle  of  virtual  veloc 

is  a  kinematically  admissible  time  independent  velocity  field.  Let  the  assoc: 

♦ 

strain  rates  be  .  Suppose  now  that  the  model  is  subjected  to  forces  which 
over  the  time  interval  0  <  t  <**,  causing  accelerations  u^  and  stresses  Q^.  For 
in  the  interval 


J 


S 


P.  ♦.dS  - 
1  1 


m  u.  ♦ .dS  = 
1  1 


iQ.  q*dS 
‘  J  J 


(36) 


Integrate  from  time  t  =  0  to  t  =  T  .  Let  the  initial  velocities  of  the  model  be  ze 

* 

and  the  final  velocities  v.  . 

1 


Js 


( 


Q.q  *  dS 

J  J 


(37) 


Advantage  has  been  taken  of  the  fact  that  is  time  independent.  We  now  suppo 


that  T-*  0  and  that 


S 


a 


o 


P.dt  -*  mv. . 
1  1 


V.  becomes  the  initial  mode  velocity  and  the  term  of  the  right  of  (37),  representing 
a  measure  of  work  done  in  the  time  interval,  will  vanish.  Finally 


1 


S 


mv.  i.dS 
1  1 


* 

mv.  ♦  dS 
1  i 


(38) 


Substituting  from  (35)  this  expression  may  be  written  either  as 


S 


♦ 

mv.v.  dS  = 
1  1 


* 


mv.  V 
1  i 


dS 


(39a) 


or 


a 


mv.  ♦  .dS 
1  1 


m  ♦.  .  dS 

1  1 


r39b) 


Thus,  referring  back  to  equation  (28),  the  velocities  for  a  mode  solution  cai 
be  written  in  terms  of  an  initial  velocity  field  and  a  time  function;  the  time  function 
[  equations  (33)  and  (34)]  and  the  initial  velocity  field  can  be  determined  completely 
in  terms  of  an  assumed  mode  shape  ♦.  and  the  physical  properties  of  the  structure. 

The  application  of  this  method  will  be  illustrated  by  the  case  of  the  beam 
shown  in  Fig.  9.  Assume  a  mode  shape  of  the  same  form  as  shown  in  Fig.  9(c)  i.e. 
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of  the  same  form  as  was  used  in  the  previous  approximation.  Thus  $  is  given 


«  S  -i|L.  0  <  X  <  i/2 


S  (x  *  0  ^  X  <  1/2 

1  o 


(40b 


The  initial  conditions  for  the  actual  problem  (Fig.  9a)  are  given  by 


V.  —  0 
1 


0  <  X  4(1  -  b)/2 


V.  -•  V 

1  o 


(1  -  b)/2  4  X  4  i/2 


To  determine  u  we  use  equation  (39b).  This  gives 
o 


mv  (— -~)  dx 
o  1 


2x  2 

m  dx 


3  vb 
^  1^ 


(21  -  b) 


Comparing  this  result  with  equation  (25),  it  is  seen  that  the  process  carried  ou 
is  equivalent  to  matching  the  initial  angular  momentum  of  half  the  beam.  The  j 
method  has  an  advantage,  however,  in  that  it  may  be  easily  applied  in  situation 


angular  momentum  cannot  be  clearly  understood.  The  time  t^  which  elapses 


the  model  comes  to  rest  is  given  by  substitution  into  equations  (3Z)  and  (34) 


‘f  = 


,  2x  .2  «  2 

m  (u^  )  dx 


4Mo 


X 

mi  V 


-2-  [1  -(1 


The  initial  velocity  at  the  center  of  the  beam  is  u^  ,  and  this  velocity  changes 
linearly  with  time.  Hence  the  central  displacement  of  the  model  is  given  by 


.  *  * 

*  o  f 


mb^  V  ^  - 

J - SL  [-i-(2 - 

12  M  ^  16  '  i  '  ^ 
o 


Comparing  this  result  with  equation  (26)  which  gives  the  displacement  obtained  from 
tlie  earlier  model,  we  find  that  the  second  approximation  is  twice  as  large  as  the 
first.  This  is  a  substantial  difference.  Without  consulting  the  correct  answer 


given  in 


^2  of 


the  paper,  we  would  expect  the  second  result  to  be  more  reliable 


than  the  first,  since  i'  involves  fewer  decisions  on  the  behavior  of  the  model. 
However,  this  is  not  necessarily  true.  Further,  the  choice  of  a  mode  shape  itself 


is  an  arbitrary  decision.  Alternative  mode  shapes  for  the  problem  under  discussio 
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are  shown  in  Fig.  10.  It  may  be  possible  that  these  shapes  would  give  better 
approximations  than  Fig.  9(c),  at  least  for  certain  values  of  b/i. 

If  the  technique  used  to  obtain  the  second  approximation  is  to  be  useful  in 
problems  where  the  exact  answer  is  not  known,  some  method  must  be  found  which 
can  differentiate  between  a  good  choice  of  mode  shape  and  a  poor  choice.  In  the 
following  sections  we  shall  attempt  to  show  how  the  ideas  used  in  establishing 
uniqueness  for  the  problems  under  discussion  can  be  used  to  make  such  a  differenti 


4.  CONVERGING  SOLUTICISS 

Using  an  argument  closely  related  to  that  used  to  prove  uniqueness,  a  rela¬ 
tion  between  the  solutions  to  two  independent  impulsive  loadings  (or  initial  velocity 
distributions)  on  the  same  structure  can  be  written,  provided  that  the  boundary 
conditions  are  the  same. 

We  consider  two  identical  structures  with  identical  boundary  conditions  (i.e, 

velocities  prescribed  zero  on  S^.  tractions  prescribed  zero  on  S,^,).  Let  the  initial 

♦ 

velocities  be  respectively  v^  (S)  and  v^  (S),  and  let  the  solutions  be  respectively 

^  ^  ^  ^ 

u.  (S,  t),  Vi.,  q.,  Q.  and  u,  (S,  t),  u.  ,  q.  ,  Q  .  As  in  equation  (4),  we  may  write, 

1  1  j  1  '  1  J  J 

A  r 

-  \  m  (u.  -  *a.  )  (u,  -  ii.  )  dS  =  \  (Q.  -  Q.  1(4  “  4  )  £iS>  0  (4! 

Jo  1  1  '  r  1  J  e  J  J  J  J  ~ 


(46a) 


where 


(u.  -  u.  )  (u.  -  u.  )  dS 

2  1  1  1  I 


(46b) 


The  initial  value  of  A  ,  at  time  t  =  0,  is  given  by 


m  * 

-rr-  (v,  -  V.  )  (v.  -  v.  )  dS 
2  1  11  1 


It  is  clear  from  (46b)  tlia*^  A  is  a  non-negative  scalar  function  of  time,  and,  from 
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(46a).  A  decreases  with  increasing  time.  The  rate  of  decrease  of  A  is  specified 

❖ 

the  right  hand  side  of  equation  (45),  and  can  be  zero  only  when  either  Q.  =  Q.  or 
•  ^ 

q.  =  q.  at  each  point  in  the  structure.  Furthermore  A  is  a  measure  of  the  differ< 
between  the  two  velocity  fields  and  ,  and  is  zero  only  when  the  two  fields  an 
identical.  According  to  this  measure,  therefore,  the  two  solutions  approach  each 
other  as  time  progresses. 

The  convergence  of  the  two  solutions  suggests  a  means  of  approximating 

unknown  solution.  If  tlie  response  to  the  velocities  v.  are  not  known,  but  the  resp 

to  velocities  v.  are  known,  we  are  assured  that  the  difference  between  the  respor 
1 

(as  measured  by  A  )  will  decrease  with  time.  If  the  difference  is  initially  not  larj 
i.e.  if  is  not  large,  we  may  expect  that  the  solutions  will  not  be  greatly  diffen 
particularly  away  from  the  initial  instants.  Furthermore,  we  would  not  expect  gr 
differences  in  the  velocities  at  any  particular  point  on  the  structure  (as  opposed  tc 
integrated  effect  represented  by  A). 

These  ideas  may  be  applied  to  attempts  to  approximate  responses  by  mea 

of  mode  solutions,  as  discussed  in  the  previous  section.  However,  in  order  to 

.  * 

substitute  a  mode  solution  for  u^  in  equation  (45)  it  must  be  a  complete  solution  tc 
an  initial  value  problem.  The  mode  solutions  used  in  the  previous  section  were 
required  to  be  kinematically  admissible  only.  If  it  is  possible  to  associate  with  th 
mode  solution  a  statically  admissible  stress  field  which  everywhere  satisfies  the 
yield  condition,  then  the  mode  solution  is  a  full  solution  to  an  initial  value  problen' 


and  may  be  substituted  into  equation  (45). 

If  the  mode  solution  is  given  as  in  equation  (28), 


(48) 


ii.*  (S.  t)  =  v.*  (S)  T  (t) 


the  accelerations  are  given  by 


*  *  dT 

u.  (S,  t)  =  v.  (S)-^(l) 


dX/dt  may  be  found  from  equation  (32).  Further,  if 


V.  =  a  4’. 
1  1 


as  in  equation  (35),  where  ♦  ^  is  a  mode  shape  and  O  is  an  arbitrary  magnitude,  if 

♦ 

may  be  shown  by  substitution  into  equations  (32)  and  (49)  that  u.  (S,  t)  is  independ 

of  O.  Tliis  means  essentially  that  in  the  mode  model  the  accelerations  are  indepei 

of  the  initial  velocity,  i.e.  the  initial  amplitude  assigned  to  the  mode.  Using 

•  * 

d'Alembert's  princple  and  considering  the  inertia  forces  (-mu.  ),  we  have  a  static 

problem  with  tractions  (in  this  case  inertia  forces)  given  on  S.^,,  and  displacementf 

20 

given  zero  on  S^.  Using  the  limit  theorems  of  Drucker,  Greenberg  and  Prager  , 
if  any  distribution  of  stresses  can  be  found  which  is  in  equilibrium  with  the  inertia 
forces  and  which  does  not  violate  the  yield  condition,  then  there  is  associated  with 
the  mode  solution  a  safe,  statically  admissible  stress  field.  Again,  this  remains 
true  for  arbitrary  a  in  equation  (50).  It  will  be  convenient  to  refer  to  a  mode  solu 
which  can  be  associated  with  a  safe  statically  admissible  field  as  a  safe  mode  solu 
If  a  solution  of  the  form  of  (48)  and  (50)  is  to  be  substituted  into  (45);  it  it 
necessary  to  determine  the  initial  mode  amplitude  a.  Since  t  [equation  (46)] 


D.  C.  Drucker,  H.  J.  Greenberg  and  W,  Prager,  "Extended  limit  design  th.eorc 
for  continuous  media,"  Quart.  App.  Math.  9,  p.  381,  1952. 
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measui’es  the  difference  between  the  solutions,  it  is  logical  to  choose  a  such  that 
A  °  is  as  small  as  possible.  Tliis  may  be  written  as 


(v.  -  a  4* .)  (v.  -  o  4.)  dS 

M  1  11  A 


Differentiating  and  solving  for  a,  we  obtain 


m  V  4  dS 
1  i 


m4>  4  dS 
1  i 


Comparison  with  equation  (39)  shows  that  the  result  is  identical  to  that  obtained 
earlier  by  a  different  argument.  Multiplying  through  by  a,  we  see  that,  using  (50), 


I  «  ❖  1  * 

mo4,a4.dS=  \  mv.  v  dS  =  \  mv.v.  dS 

S  ^  ^  J  S  ^  ^  4  ^  ^ 


9^  ^ 

m  (v.  -  V.  )  V.  dS  =  0 
111 


Equation  (53a)  may  be  used  to  show  that,  if  (5Z)  is  satisfied, 


=  \  v.v.  dS 

Os  ^  ‘  ‘ 


m  *  * 

V.  V.  dS 
2  11 


i.e.  the  initial  measure  of  difference  is  simply  the  initial  difference  between  the 

energies  of  the  solutions.  Equation  (53b)  shows  that  the  difference  between  and 

is  orthogonal  to  ;  this  result  is  suggestive  of  elastic  normal  mode  analysis, 

* 

considering  to  be  analogous  to  any  one  m.ode. 

This  discussion  indicates  that  the  choice  of  a  safe  mode  solution  to 

approximate  the  unknown  solution  leads  to  initial  conditions  identical  to  those  found 

earlier,  and  further  to  an  assurance  that  a  measure  of  the  differences  between  the 

solutions  must  decrease.  If  the  initial  difference  is  small,  i.e.  small  compared 

to  the  initial  energy  in  any  one  solution,  we  are  provided  with  a  reliable  check  of 

the  validity  of  the  one  degree  of  freedom  model.  No  arbitrary  assumptions,  other 

than  the  choice  of  a  mode  shape,  need  be  made. 

One  further  useful  piece  of  information  may  be  added.  From  a  proposition 
21 

due  to  Martin  ,  the  time  t  at  which  the  velocities  vanish  in  the  real  solution 
can  be  written 


(55) 


▼  * 

where  v.  is  the  initial  mode  velocity  and  D  is  the  initial  rate  of  energy  dissipation 
in  the  mode.  (This  proposition  does  not  require  that  the  mode  be  a  safe  mode.) 
However,  from  (34),  (32)  and  (53a), 


J.  B.  Martin,  "Impulsive  loading  theorems  for  rigid  plastic  continua,  "  Proc. 
Eng.  Mech.  Div. ,  A.  S.  C. E .  ,  90  (EM5),  p.27,  1964. 


mv.  V.  dS 
i  i 

1 , 


mv.v.  dS 
1  1 


D  ds 


(56) 


provided  that  the  initial  amplitude  of  the  mode  satisfies  the  optimum  requirement 
that  be  a  minimum.  Hence 


(57) 


i.e.  the  velocities  in  the  model  vanish  before  the  velocities  in  the  actual  solution. 

* 

Experience  shows  that  t^  and  t^  are  very  often  equal,  and  in  most  cases  the 
difference  is  small. 

In  tlie  following  section  we  shall  discuss  the  application  of  this  method  to 


illustrative  examples. 


5.  EXAMPLES 


In  order  to  demonstrate  the  effectiveness  of  the  approximation  using  the  mode 
solution  we  shall  compare  the  results  of  the  solutions  discussed  in  ^ 2  with  their 
approximations . 


(i)  Consider  first  the  clamped  beam  shown  in  Figs.  6(a)  and  9(a).  The  final 

central  displacement  is  given  in  equations  (17)  and  (18),  and  is  a  function  of  b/i. 

A  mode  solution  has  been  found  for  this  problem,  Fig.  9(c),  and  the  final  displacement 
is  given  in  equation  (44).  It  remains  to  check  that  the  mode  chosen  is  indeed  a  safe 
mode,  and  to  find  the  relative  value  of  the  initial  difference 

The  acceleration  of  the  central  point  in  the  beam  in  the  mode  solution  may 
be  easily  calculated  from  the  initial  central  velocity  [  eqja.tion  (42)]  and  the  time  of 
duration  [  equation  (43)]  since  it  is  known  that  the  velocity-time  relation  is  linear. 

This  is  given  by 


..  * 

u 


24M 

o 


(5=*) 


Using  d'Alembert's  principle,  the  loading  diagram  for  the  beam  is  shown  in  Fig.  11(a) 
and  the  bending  moment  diagram  is  shown  in  Fig.  11(b).  It  is  readily  shown  that  the 
maximum  bending  moment  is  M^,  and  that  it  always  occurs  in  the  center  of  the  beam. 
Hence  the  mode  is  a  safe  mode  solution.  In  order  to  compute  A®  we  make  use  of 

equation  (54).  A  relative  measure  of  the  magnitude  of  may  be  obtained  by  dividing 

o  o 

&  by  the  initial  energy  of  the  actual  problem.  Let  this  be  K  .  Then 
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o  1  2 

K  s  — r-  mbv 
2  o 


mbv 


2 


dx 


mbv 


miv 


2b 


-  (-f  )^]  ^ 


,  (59a) 


(59b) 


(59c) 


The  final  central  displacement  computations  have  been  plotted  in  Fig.  12.  In  ad< 

to  the  actual  and  mode  solutions,  the  upper  bound  which  may  be  computed  by  a  m 

21 

proposed  by  Martin  is  given.  In  Fig.  13  the  mode  solution  error  and  the  ratio 
A  are  plotted  on  the  same  figure.  This  diagram  shows  that  the  error  is  sm 

when  A°/K°  is  small,  and  that  the  error  increases  as  A°/K°  increases.  It  cai 
be  expected  that  the  error  is  zero  when  A°/K°  is  a  minimum,  and  indeed  this  is 
so. 

Figs.  14  and  15  give  further  details  for  the  case  of  h/i  -  1.  Fig.  14  gi 
the  central  velocity-time  curves  for  the  actual  solution  and  the  mode  solution.  It 
be  seen  that  the  velocities  coincide  after  one  third  of  the  total  deformation  time  h 
elapsed.  In  fact  the  velocities  coincide  everywhere  on  the  beam  during  this  peric 
and  L  [  equation  (46b)]  is  zero.  Thus  t  decreases  rapidly  (and  in  this  case 
vanishes)  even  though  (Fig.  13)  has  an  appreciable  magnitude. 


Fig.  15  compares  the  final  shapes  of  the  deformed  beams  for  the  case 

b/f  s  1,  showing  that  the  displacements  in  the  mode  solution  overestimate  the  actual 

displacements  in  some  regions,  and  underestimate  them  in  others. 

A  mode  solution  for  the  pinned  end  case  of  the  beam  discussed  above 

(Fig.  7a)  may  be  obtained  in  a  manner  identical  to  that  for  the  clamped  case.  The 

mode  shape  is  taken  to  be  that  in  Fig.  9(c).  The  initial  amplitude  is  not  affected  by 

the  change  in  boundary  conditions,  but  the  deformation  time  t^  is  simply  doubled. 

It  may  readily  be  shown  that  this  is  a  safe  mode  solution.  A  comparison  of  the 

19 

approximate  solution  with  the  actual  solution  for  the  central  displacement  for 
two  values  of  b/f  is  given  in  Fig.  16. 

(ii)  Consider  secondly  the  cantilever  with  an  attached  tip  mass  shown  in  Fig.  4. 

The  solution,  giving  the  tip  mass  velocity  as  a  function  of  time,  appears  in  equations 
(8)  and  (10).  The  final  tip  mass  displacement  is  given  in  equation  (11). 

In  this  case  a  good  guess  at  a  mode  solution  is  clearly  a  velocity  field 
involving  a  rigid  body  rotation  about  the  support,  shown  in  Fig.  17.  The  velocity 
field  and  rotation  rate  at  the  base  may  be  written  in  terms  of  z  ,  the  velocity  of 
tlie  tip  mass.  From  equations  (28)  and  (32),  the  acceleration  of  the  tip  mass  will  be 


..  ♦ 

z 


z  (M^z  /i) 


M 


.  *  2 

G(z  )  + 


m  (“^  z*)  ^  dx 


Gf  (1  +— Y) 


(60) 


''o 

where  y  -  ml/ZG  as  before.  The  bending  moment  diagram  may  now  be  drawn,  and 
is  shown  in  Fig.  18.  It  may  readily  be  seen  that  the  bending  moment  does  not  exceed 


at  any  point  for  any  value  of  Y  •  Hence  the  chosen  approximation  is  a  saf< 
mode  solution. 

*  * 

The  optimum  initial  value  of  the  tip  velocity  z  may  be  found  as  befor 

* 

The  actual  initial  velocity  is  zero  except  at  the  tip  where  it  has  value  v.  If  v 

.  ♦ 

is  the  initial  value  of  z  ,  we  require,  from  equation  (53a), 


i 

X  ^2. 
m  (-^  V  )  dx 

o 


The  tip  mass  velocity  is  then  given  by 


where,  from  equation  (56) 


M  ( 
o 


G/v 

M 

o 


(62) 


Fig.  19  shows  the  tip  mass  velocities  plotted  for  the  actual  solution 
[  equations  (8)  and  (10)]  and  for  the  mode  solution,  for  the  particular  case  Y  =  1. 

In  this  case  it  is  seen  that  the  velocities  are  the  same  for  M  t/Giv  >0.  167:  it  ma 

o  — 

readily  be  shown  that  A  vanishes  at  time  M^t/Gfv  =  0.  167. 
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The  final  tip  displacement  is  a  function  of  the  parameter  y  ,  and  hence  the 
accuracy  of  the  approximation  will  depend  on  Y  .  The  final  tip  displacements  have 
been  plotted  as  a  function  of  Y  in  Fig.  20,  and  show  that  the  discrepancy  varies 
from  zero  at  Y  =  0  to  an  underestimate  of  about  15%  at  y  =  3.  On  the  same 
figure  has  been  plotted.  When  y  =  3  /K*^  has  the  value  0.  67.  Despite 

this  large  value  of  the  approximation  is  fairly  good. 

One  further  example  is  summarized  in  Fig.  21.  In  this  case  the  energy 

.0/0 

ratio  ^  /K  is  extremely  small.  Considering  the  areas  under  the  curves  to 
obtain  displacements  at  the  center  of  the  beams,  it  can  be  seen  that  the  difference 
between  the  actual  solution  and  the  mode  solution  is  extremely  small. 


6.  CONCLUSIONS 


In  this  paper  an  attempt  has  been  made  to  rationalize  the  setting  up  of  c 

y 

one  degree  of  freedom  approximation  in  elementary  rigid-plastic  theory  for  imp 
loading.  The  method  requires  that  a  mode  shape  be  chosen;  thereafter,  the 
deceleration  (or  the  equivalent  spring  force  and  mass)  and  the  initial  mode  veloc 
follow  without  further  assumptions.  The  concept  of  a  safe  mode  and  the  analysii 
given  in  Section  4  provide  a  criterion  by  which  good  mode  approximations  may  b 
recognized.  For  illustrative  purposes  the  method  has  been  applied  to  extremely 
simple  examples.  There  is  no  conceptual  difficr.lty  in  applying  the  technique  to  1 
more  complicated  cases. 

It  is  seen  that  even  in  cases  where  is  fairly  large  the  approximj 

tion  can  be  reasonably  good.  Only  where  jf/K^  is  small  however,  can  it  be 
taken  that  the  approximation  will  certainly  be  good.  The  other  cases  emphasize 
that  good  approximations  can  be  obtained  whe;  is  large,  as  undoubtedly 

good  approximation  could  be  obtained  with  mode  solutions  wjiic!;  were  .not  safe  ir. 
the  sense  used  in  this  paper.  This  paper  shows  o  ily  that  a  certain  limited  class 
of  mode  approximations  can  be  considered  reliable. 

The  techniqi..e  described  could  be  applied  with  c'  anges,  to  problems  wi1 
time  dependent  loading  and  to  certain  oti  er  viscous  type  materia!  idealizations 
such  as  rigid -visco-plastic .  Attempts  are  being  made  to  develop  u.sefiil  approxi¬ 
mating  techniques  for  these  cases. 
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II.  AMTAACT 


An  attempt  is  made  to  provide  a  rational  method  of  constructing 
one  degree  of  freedom  approximations  for  impulsively  loaded  metal 
structures  which  are  analysed  according  to  an  elementary  rigid-plastic 
theory.  The  approximation  follows  automatically  from  a  chosen  mode 
shape,  and  a  criterion  for  determining  good  mode  shapes  is  introduced. 
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